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We have studied the magnetoresistance in a high-mobility Si inversion layer down to low electron 
concentrations at which the longitudinal resistivity pxx has an activated temperature dependence. 
The angle of the magnetic field was controlled so as to study the orbital effect proportional to the 
perpendicular component B±^ for various total strengths Btot- A dip in p^x, which corresponds 
to the Landau level filling factor oi v — A, survives even for high resistivity of pxx ~ 10** Vl at 
T = 150 mK. The linear Btot-dependence of the value of B± at the dip for low Btot indicates that 
a ferromagnetic instability does not occur even in the far insulating regime. 

PACS numbers: 73.40.Qv, 71.30.+h, 73.20.Qt 
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I. INTRODUCTION 



There has been great attention to the fundamental 
properties of strongly correlated two-dimensional (2D) 
electron or hole systems in the last decade, due in part to 
the discovery of the zero-magnetic-field metal-insulator 
transition (MIT)i2i^ The strength of the Coulomb in- 
teraction between the carriers is characterized by the 
Wigner-Seitz radius r^, which is equal to the ratio 
of the Coulomb energy per electron to the Fermi en- 
ergy. In the Fermi liquid theory, the effective mass 
m*, the effective factor g* and the spin suscepti- 
bility X* oc g*m* are renormalized by r^, and they 
are expected to be enhanced as increasesA^ The 
dimensionless parameter can be written as Ts — 
TT^^'^{e/h)'^{mb/Keo)Ns~^^^ and controlled by changing 
the electron concentration Ng, where nii, is the band 
mass and k is the average dielectric constant. The 
enhancement of the spin susceptibility with increas- 
ing Ts (decreasing Ng) has been observed in various 
2D electron systems (2DESs) formed in Si metal-oxide- 
semiconductor field-effect-transistors (M0SFETs)p£i2*a 
GaAs/AlGaAs,i" Si/SiGe" and AlAs/ AlGaAsi^ het- 
erostructures. For 2DESs without disorder, the ferro- 
magnetic transition is expected to occur at rg ^ 26 be- 
fore the Wigner crystallization at r^. 35}^ Recently, 
the divergence of x* at or near the MIT with Ts ~ 9 has 
been reported for Si-MOSFETsji^ while the results in 
Refs. andlTil do not support the occurrence of a ferro- 
magnetic instability at the MIT. Another 2D Fermi liquid 
system, '^He absorbed on graphite, also shows a tendency 
for X* to diverge as the system goes into localizationiiLi^ 

A possible ground state in the insulating regime of 
high-mobility Si-MOSFETs is a pinned Wigner crystal 
(WC) or glass. Pudalov et al. observed nonlinear dc 
conduction with a sharp threshold electric field in the in- 
sulating regime of Si-MOSFETs and attributed it to that 
of a pinned WC; ^- while it was also discussed in terms of 
the single particle localization picture taking into account 
a Coulomb gapi22i2i Chui and Tanatar found from their 
Monte Carlo studies that the WC can be stabilized at r. 



as low as 7.5 in the presence of a very small amount of 
disorder in the oxide layer i22iS^ In the WC at rather low 
Tg, exchanges among neighboring electrons are expected 
to occur frequently. The amplitudes of several types of 
ring exchanges in Si-MOSFETs were calculated^'i'-?^ to 
be in the order of 0.1 K at ~ 8 using the WKB 
method developed by Rogeri^ The strength of the fer- 
romagnetic interaction, which arises from exchanges of 
an odd number of particles, is comparable to that of the 
antiferromagnetic interaction from exchanges of an even 
number of particles.-^- Furthermore, the valley degree of 
freedom in Si inversion layers makes the system more 
complicatediSl It seems hard to predict theoretically the 
magnetic ground state of the WC in Si 2DESs. 

In this work, we have performed systematic mea- 
surements of magnetoresistance of a high-mobility Si- 
MOSFET in order to study the electronic and spin states 
in the insulating regime. The perpendicular component 
B± of the magnetic field was controlled by rotating the 
sample for various total strengths Btot so as to investigate 
the orbital effect independently of the Zeeman effect. We 
observed that a dip in the longitudinal resistivity p^ji, 
which corresponds to the Landau level filling factor of 

— 4:, remains even in the far insulating regime where 
the Landau levels are expected to be smeared out. The 
value of B± at the dip shows a linear iJtot-dependence 
in the low Btot region below a kink indicating the onset 
of the full spin polarization. The results strongly sug- 
gest that a ferromagnetic instability does not occur in 
the insulating regime. 



II. SAMPLE AND EXPERIMENTAL METHOD 

We used a (OOl)-oriented Si-MOSFET sample with 
a peak electron mobility /ipoak = 2.4 m^/Vs at Ng = 
4 X 10^^ m~^ and T — 0.3 K. It has a Hall-bar geometry 
of total length 3 mm and width 0.3 mm. The estimated 
Si02 layer thickness is 98 nm. Standard DC four-probe 
techniques were used to measure pxx and the Hall resis- 
tivity pxy The potential probes were separated by 1.5 
mm and the excitation voltage had been kept 0.4 mV to 
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FIG. 1; Longitudinal resistivity as a function of B± at Btot = 
9 T and T = 150 mK for various Ns indicated in units of 
The inset shows the data at T = 180 mK for 
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ensure that measurements were taken in the I-V hnear 
regime. The electron concentration Ns was controlled by 
varying the gate voltage and determined from the Hall 
coefficient measured at T = 3 K. The MIT is observed at 
the critical electron concentration Nc = 0.97 x 10^^ m~^ 
in the absence of the magnetic field. It is estimated that 
rs = 8.4 at Ns = Nc with mb = 0.19mo and k = 7.7^ 
The sample was mounted on a rotatory thermal stage in a 
dilution refrigerator together with a GaAs Hall generator 
and a Ru02 resistance thermometer calibrated in mag- 
netic fields. The rotatory thermal stage was cooled via a 
silver foil linked to the mixing chamber and the temper- 
ature was accurately controlled by a heater on the stage. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 



Figure 1 shows the B^-dependence of the longitudi- 
nal resistivity pxx at Btot = 9 T for various Ns with 
a wide range of pxx- As discussed later, the electron 
spins are expected to be fully polarized at Btot = 9 T 
for these values of A''. For TV. = 1.97 x 10 
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ima in p^x are observed at the Landau level filling factors 
1/ ~ 4, 6, 8 and 10, which correspond to the Shubnikov- 
de Haas (SdH) oscillations of the spin polarized system. 
Note that twofold valley degeneracy remains in Si(OOl) 
inversion layers. As Ns decreases, pxx drastically in- 
creases and the dips at « 6, 8 and 10 are smeared 
out. It is reasonable that the SdH oscillations disap- 
pear when a dimensionless parameter lOcT becomes less 
than unity. Here u>c is the cyclotron frequency and r 
is the cyclotron scattering time. This condition can be 
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FIG. 2; Data in the insulating regime for Bt, 



6 T and 



Ns = 1.02 X 10^^ m~^. The electron spins are expected to be 
fully polarizedi- (a) Bx-dependence of p^x for various tem- 
peratures, (b) Arrhenius plots of pxx for B± = T, 0.94 T 
and 1.19 T (at the dip), (c) Activation energy as a function 
of B±. 



rewritten as p|| > h/ve^ (« pxy) if t is replaced by the 
classical scattering time obtained from the resistiv- 
ity p|| = Pxx{B±_ = 0). Measurements on low-mobility 
Si-MOSFET&29 and GaAs/AlGaAs 2DES92a have shown 
that T is comparable or smaller than Tc. Thus, in a simple 
picture, the SdH oscillations are not expected to appear 
for p|| 10** 57. However, the dip at ^' « 4 remains for 
very high resistivity up to pxx ~ 10^ ^2. Distinct pxx 
minima &t v k, integer in the insulating regime of high- 
mobility Si-MOSFETs have also been observed for v \ 
and 2 in perpendicular magnetic fields {B±_ — Btot)^^^ 
It is well known that the usual SdH oscillations origi- 
nate from the i?_L-dependence of the density of states 
at the Fermi level ep- The longitudinal conductivity 
cTxx = Pxx/iPxx'^ + Pxy^) has a minimum when sp lies 
in a gap between Landau levels. For pxx ^ Pxyi a mini- 
mum in a XX leads to a minimum in pxx- In the insulating 
regime with pxx 3> Pxyi on the other hand, it leads to a 
maximum in pxx^ which is contrary to the experimental 
results for v w integer. 

Figure 2(a) shows the i?_L-dependence of pxx in the in- 
sulating regime for various temperatures. Pxx decreases 
drastically with increasing temperature even for the min- 
imum sX u K, A while pxx increases with T for minima in 
the usual SdH oscillationsiS While the dip at « 4 in 
Fig. 2(a) is gradually smeared out as T increases, the po- 
sition of the dip does not depend on T. Figure 2(b) shows 
Arrhenius plots of pxx for different B±^ . The dashed lines 
are least-square fits to the experimental data and repre- 
sent Pxx = Po&xp{Ea/2T).^^ The i3j_-dependence of pxx 
at low temperature is attributed to a change in the ac- 
tivation energy Ea rather than that in the prefactor pQ. 
The obtained Ea is shown in Fig. 2(c) as a function of 
B±. In the low B± region, Ea decreases almost linearly 
with increasing B±. This might be associated with the 
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delocalization effect of the magnetic field in the strongly 
localized regime. Ej^ shows a dip at « 4. Although 
it might be a trace of the Landau level formation, the 
depth of the dip of AEa ~ 0.1 K is much smaller than 
the Landau level spacing of TilUc ~ 8.5 K at B± = 1.2 T 
expected from the band mass of mf, = 0.19mo. At this 
stage, the origin of p^x or Ea minima at w integer 
observed in the insulating regime is not understood. We 
have investigated i?_L-dependence of pxx also for other 
2D systems with high resistivity of Pxx ^ 10^ $7 at low 
temperatures down to 100 mKf24 Dips at « 1 and 2 are 
observed for GaAs hole systems with rg ^ 10, while an 
insulating GaAs electron system with small (« 3) only 
shows a broad minimum resulting from negative magne- 
toresistance at low B± and positive one at high B± owing 
to the shrinkage of the electron wave function. Such B±- 
dependence of Pxx in a GaAs 2DES was also observed 
in Fig. 2 of Ref. The collective motion of electrons 
might cause the dips in pxx Bd, v ~ integer observed in the 
insulating regime of the strongly correlated 2D systems. 

The values of B± at pxx minima depend on i?tot- 
Typical data in the low resistivity region are shown in 
Fig. 3(a). The positions of the pxx minima at w 4 and 
6 in the SdH oscillations shift toward low-i?^ side as Stot 
decreases. This can be explained as the result of a de- 
crease in the fraction of "spin-up" elect ronsAiiSLSLSS 
Similar behavior is observed in the insulating regime. 
In Fig. 3(b) and (c), the data for Ns = 1.02 and 
0.94 X 10^^ m~^ are shown, respectively. The dip at « 4 
shifts toward low-i?j^ side as Btot decreases. In Fig. 4, 
the positions of the pxx minima determined taking into 
account the negatively i?j^-dependent baseline are shown 
as l/t'min — eB±/hNs for different Ng. Overall behav- 
ior for low Ns is similar to that for high Ns, i.e., l/i'min 
increases almost linearly with _Btot before i3tot exceeds a 
critical value Be indicated by arrows. We consider that 
Be corresponds to the onset of the full spin polarization 
of 2D electrons. The dotted lines represent = 4 or 
6 assuming that the spin polarization P ~ 2Ni^/Ns — 1 
increases linearly with Btot for -Btot < Be- Here and 
iV-f are the Landau level filling factor and the concentra- 
tion of spin-up electrons, respectively (i^| — hN^/eB±). 
The SdH oscillations depending on are also observed 
in a Si/SiGe sample with much lower resistivity. ^"'^ Al- 
though the reason for the survival of the dip in the far 
insulating regime is unclear, the dotted line for = 4 
well reproduces the experimental results for Btot < Be- 
The i?tot-dependent behavior for Btot < Be strongly sug- 
gests that the spin polarization is not completed and a 
ferromagnetic instability does not occur. 

A gradual increase in 1 / ^'min with i?tot is observed 
above Be for which the spin polarization is expected to 
be completed. Similar behavior is found for pxx min- 
ima at « 1 and 2 in the insulating regime as shown 
in Fig. 5. We estimate that the change in Ns due to 
the B-dependent shift of the chemical potential of the 
2DES for a fixed gate voltage is negligible (in the order 
of 10~^Ns). It seems possible that spin-down electrons 
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FIG. 3: _Bj_-dependence of pxx for various Stot. (a) Data in 
the low resistivity region with Ns — 1.51 x 10^^ m~^. (b)(c) 
Data in the higher resistivity region with A'^s — 1.02 and 0.94 x 
10^^ m~^, respectively. 



occupy deep levels due to impurity potentials even at 
Btot = Be and are gradually released into non-trapped 
states with an up-spin at higher magnetic fields. In the 
Si/SiGe sample,-^- the increase in l/t'min with Btot is not 
observed above Be- 

The critical magnetic field Be is also obtained from a 
kink in a magnetoresistance curve in the in-plane mag- 
netic field B\\ since the By-dependence of pxx is asso- 
ciated with the spin polarizationiSiiiiSSiSLSS Figure 6(a) 
shows the data in the insulating regime. While pxx has 
strong temperature dependence, the -Bii -dependence at a 
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FIG. 4: Position of pxx minima 1/z^min = eB±/hNs. Data 
for (a) Ns = 0.94 x 10^^ m~^ (b) 1.02 x 10^^ m-^ (c) 1.27 x 
lO^'^ m-2 and (d) 1.51 x lO^'^ m ^ are shown. The arrows 
indicate Be- The dotted lines represent i^^ = 4 or 6 (see text). 



constant temperature shows a steep increase in the low 
region and a saturation in the high region. The 
critical magnetic field determined from the magnetoresis- 
tance curve is consistent with that obtained from Fig. 4 
as shown later in Fig. 7. It also suggests that a ferromag- 
netic instability does not occur in the insulating regime. 
In Refs. IT3 and lTsl the authors claimed that Be tends to 
vanish at an electron concentration close to the MIT at 
B = based on the analysis of the -dependence of pxx 
obtained in the metallic side. However, the magnetore- 
sistance in the far insulating regime was not studied. 

In Fig. 6(b), we show the activation energy Ea de- 
termined from the Arrhenius temperature dependence of 



1 — I — I — I — I — I — I — I — I — 
7 = 150 mK 
/V,=0.86x10'=m ' 




2.5 3.0 

e,(T) 



4.0 



FIG. 5: Bx-dependence of p^x at T = 150 mK and Ns = 
0.86 X lO^^m"'^ for various -Btot. pxx minima at ~ 1 and 2 
are observed. The inset shows the position of the pxx minima. 



10' 



^10%^* ■ 



j/^* iAiAiAii AAAAAAAAAA 



(a) 




T = 300 mK 



• 0.82 

* 0.86 
■ 0.94 



(b)- 



aE^/ae|| = 2.9i|iB 
aE^/ae|| = 2.37no 



• 0.82 
■ 0.94 



6 



10 



eii(T) 



FIG. 6: (a) In-plane magnetic field dependence of pxx for dif- 
ferent A^s at T = 300 mK. The arrows indicate Be determined 
from Pxx vs B|| data obtained for various temperatures, (b) 
Activation energy determined from T-dependence of pxx for 
Ns = 0.82 and 0.94 x 10^^ m'^ 



Pxx ■ If we assume Ea as an energy gap for an elementary 
excitation, the magnetization change due to the excita- 
tion can be obtained thermodynamically via the relation 
SM ~ —dEA/ dB\\. A similar method was used for the 
study of the energy gap for the odd-integer quantized 
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FIG. 7: Be obtained from the Btot-dependence of 1/z^min 
(closed diamonds) and that from the By -dependence of pxx 
(open circles) are shown as a ratio to the non-interacting value 
Bq. The dashed curve is an estimation from the spin suscep- 
tibility obtained by Pudalov et al^ in the limit of small mag- 
netic fields. The dotted curve represents the values of B\\/ Bq 
for a tentative boundary in the B\\-Ns plane determined from 
pxx at T = 150 mK. pc = 60 kfi is the critical resistivity at 
= Nc and B = 0. 



Hall statesi^ In our Si-MOSFET, the average distance 
of electrons from the Si/Si02 interface^^ is by one or- 
der of magnitude smaller (« 3.8 nm) than the magnetic 
length Iq = {h/eB^^y/^ in the low region where the 
steep increase in Ea is observed. Since does not cou- 
ple to the orbital motion of electrons in this case, SM 
should be attributed to electron spin flips. The slope 
dEA/dB\\ in the low B\\ region is somewhat larger than 
+2iJ,B expected from a single spin flip. Here, /is is the 
Bohr magneton and the bare (7-factor in silicon is close 
to two. The observed large \SM\ cannot be explained by 
a single particle picture. 

Figure 7 shows Be as a function of Ng. They 
are normalized by the non-interacting value Bq = 
TTh^Ns/2^B'mb with mt, ~ 0.19mo. If P increases lin- 



early with i?tot below Be, Bc/Bq is equal to the inverse 
of the ratio of x* to non- interacting susceptibility xo- 
The dashed curve represents Xo/x* obtained from the 
data by Pudalov et al^ The dotted curve represents the 
values of B\\/Bq for a tentative boundary in the B\\-Ns 
plane on which p^x aX T — 150 mK is equal to the criti- 
cal resistivity = 60 kf7 for the MIT at B = OiS While 
the magnetic-field-induced MIT occurs in the interme- 
diate Ng range, Bc/Bq is in good agreement with Xo/x* 
obtained in the metallic regime in the limit of small mag- 
netic fields.^ Bc/Bq gradually decreases with decreasing 
Ng (increasing Vg) in the whole Ng range and the iV^- 
dependence shows no distinct anomaly. 

IV. SUMMARY 

In summary, we have studied the low-temperature 
magnetoresistance in a high-mobility Si-MOSFET sam- 
ple. Even in the far insulating regime with p^x ~ 10^ ft, 
we observed a dip in pxx at « 4 in the S^-dependence 
at Btot = 9 T. The critical magnetic field for the onset 
of the full spin polarization determined from the -Btot- 
dependence of the value of B± at the dip, which agrees 
with that from the magnetoresistance curve in the in- 
plane magnetic field, indicates that a ferromagnetic in- 
stability does not occur in the insulating regime. 

It was found that -dependence and -B|| -dependence 
of Pxx at low temperature in the far insulating regime are 
the results of those of the activation energy Ea in the Ar- 
rhenius temperature dependence. However, the origins of 
the dip in Ea at « integer and the steep increase in Ea 
with _B|| for low are not understood, while they might 
be associated with strong electron correlations. Further 
investigations are required. 
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